Photoassociation spectroscopy of a Spin-1 Bose-Einstein condensate 
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We report on the high resolution photoassociation spectroscopy of a 87 Rb spin-1 Bose-Einstein 
condensate to the l g (P 3 / 2 ) v = 152 excited molecular states. We demonstrate the use of spin 
dependent photoassociation to experimentally identify the molecular states and their corresponding 
initial scattering channel. These identifications are in excellent agreement with the eigenvalues of a 
hyperfine-rotational Hamiltonian. Using the observed spectra we estimate the change in scattering 
length and identify photoassociation laser light frequency ranges that maximize the change in the 
spin-dependent mean-field interaction energy. 

PACS numbers: 34.50.Rk, 33.15.Pw 



In a spinor Bose-Einstein condensate [l|, H, d, H| , the 
delicate interplay between different atomic spin orienta- 
tions results in a small spin dependence of the collisional 
interaction energy [5j. Though the spin dependence is 
small relative to the total interaction energy (e.g. ~ 0.5% 
in 87 Rb), the coherence of the condensate allows obser- 
vation of novel phenomena such as coherent spin mixing 
and dynamics [q, as, m , metastable states [HI], s P m 
domains 12|, ll3j and quantum tunneling across spin do- 
mains [14 1 . The spin-dependent interaction energy arises 
from the small difference in the s-wave scattering lengths 
of the allowed angular momentum channels and manifests 
in anti- or ferromagnetic properties for a spin-1 conden- 
sate depending on the algebraic sign of the difference. 
In the two spin-1 condensates investigated to date, 23 Na 
and 87 Rb, the former is anti-ferromagnetic and the latter 
is ferromagnetic. 

Magnetic 



151 ] and optical [16( Feshbach resonances 



have been employed to dynamically change the scattering 
lengths which offers the tantalizing prospect of influenc- 
ing the spinor dynamics or even possibly changing the 
magnetic nature of the condensate. However, spinor dy- 
namics are suppressed for magnetic fields greater than a 
few 100s of mG @, E3]- The fields used in magnetic Fes- 
hbach resonance experiments are much larger than this 
and thus can not be used for altering the spinor dynam- 
ics. Photoassociation around single resonance lines also 
has limited flexibility since the atom loss rate is very 
high near the resonances where the change in scattering 
lengths is appreciable. To circumvent this problem it has 
been suggested to use the effect of multiple photoasso- 
ciation lines from a vibrational level with a rich hyper- 
fine structure [l^|. This allows parameters that depend 
on multiple scattering lengths to be enhanced between 
lines from different scattering channels while the atom 
loss rate is reduced. Specifically, variation of the spin de- 
pendent interaction strength could be optimized if there 
are adjacent molecular states that occur through different 
scattering channels. In order to identify the existence of 
such lines, collision channel selective spectroscopy needs 
to be performed which requires control over the Zeeman 
states of the colliding atoms. Note that although the 



effects of hyperfine levels of the colliding atoms on the 
photoassociation spectrum has been previously observed 
~ 19l 2ol ] , there has not been any collision channel selective 



spectroscopy. 

In this Letter, we report on experimental photoassoci- 
ation spectroscopy of a spinor condensate. By observing 
the effect of the Zeeman states of the colliding atoms on 
the photoassociation spectrum, we extract information 
about the angular momentum of the molecular states. 
We identify for the first time molecular states that are 
allowed for a specific scattering channel while forbidden 
for another. From the spectrum, we calculate the atom 
loss rate and the expected change in the mean-field in- 
teraction energy. Using the observed spectra and these 
calculations, we identify suitable photoassociation light 
frequencies that maximize change of the spin dependent 
interaction strength with modest atom loss rates. 

The experiment is performed on 87 Rb condensates cre- 
ated directly in an optical trap [l| . After the condensate 
is formed, it is exposed to a photoassociation (PA) laser. 
To measure the spectrum of the molecular excited states, 
the condensate population losses are measured for dif- 
ferent frequencies of the PA light. The PA excitation 
laser frequency is tuned near the l g (i"3/2) v = 152 state 
which has a rich hyperfine structure and a binding en- 



ergy of 24.1 cm 1 [21 



|22j below the D 2 line of 87 Rb. 
The PA laser has a focused waist of 80 /im, and its fre- 
quency is actively stabilized with an accuracy ^5 MHz 
using a transfer cavity locked to a stabilized diode laser. 
Fig.Q] shows the observed photoassociation spectrum of 
the l g (P3/2) v — 152 state taken using condensates with 
/ = l,m/ = —1) and |/ = l,m/ = 0) spin states. A 
portion of the \f = 1,771/ = —1) spectrum was observed 
in [22j]. The PA light is ramped up to 3.8 mW in 50 ms, 
after which it remains on for 100 ms. Successive data 
points are separated by 5 MHz. The weaker lines iden- 
tified by k through tt have been confirmed using more 
intensity and longer time but are shown here under the 
same conditions as the stronger lines for consistency. The 
inset shows a line too weak to be made out in the larger 
scan. In order to enhance the visibility of this line, the 
PA light power is increased to 11 mW and left on for 
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FIG. 1: Observed photoassociation spectrum of the l g (P3/2) v — 152 state for mj = —1 (upper) and m/ = (lower). These 
spectra are obtained after 150 ms of exposing a BEC to a PA light with 3.8 mW and 80 micron beam waist. The inset box 
for rrij = — 1 is obtained using 350 ms and 11 mW with averaging to enhance visibility of this weaker line. Plots are offset 
for clarity. The origin of the hyperfine-rotation Hamiltonian fit is used as the zero point for the plot. A stick spectrum with 
approximate \F,f,I,i) labels is given. 



300 ms, and it has been averaged over four scans. For 
each data point, an absorptive image of the condensate is 
taken 12 ms after the PA and trapping lasers are turned 
off. The condensate population is counted and normal- 
ized to an average value taken under the same conditions 
but with the PA light frequency detuned far from any of 
the molecular lines observed in Fig.Q] 

The l g potential is categorized as a Hund's case (c). 
The relevant good quantum numbers for this potential 
are total molecular angular momentum, F, total nuclear 
spin, /, and f2, the projection of the rotational angular 
momentum, J, on the internuclear axis. For this case, 
the projections Fg, lg of F&nd I on the internuclear axis 
are almost good quantum numbers, while J is not [23[. 

The possible values of F are determined by bosonic 
symmetry and angular momentum addition rules. For 
two identical spin-1 bosons, the allowed s-wave collision 
channels are F = 2 or 0. Addition of the angular momen- 
tum of the colliding atoms and the PA photon specifies 
the possible total angular momentum, F, of the available 
molecular states for each scattering channel. Photoasso- 
ciation through the total spin 2 scattering channel gives 
molecular F numbers of 1, 2, and 3 while the total spin 
scattering channel restricts F to 1. 

In the following, we use the experimental properties 
of the photoassociation spectrum to identify these quan- 
tum numbers for the observed molecular lines. The most 
striking difference between the two PA spectra in Fig.Q] 
is that the lines e and 8 appear only for the condensate 
containing the m/ = spin state. Since two mf = — 1 
atoms can only scatter through the total spin 2 channel, 
while two mf = atoms access both the total spin 2 and 



channels, this observation indicates that these lines oc- 
cur through the total spin scattering channel, restrict- 
ing F for the molecular state represented by these lines 
to 1. This conclusion is further confirmed by repeating 
the same experiment with condensates containing both 
m/ = — 1 and m/ = 1 spin states and also with conden- 
sates containing m/ = — 1 and m/ = (see Fig. [2]). These 
spin state combinations are prepared using a microwave 
manipulation technique similar to the one used in Ref. 
[7|. The populations of different mf states are counted 
by spatially separating them using a Stern-Gerlach field 
during the time of flight. An mt = — l(rnt = 1) can only 
access the total spin scattering channel through colli- 
sion with an mt — l(mf — —1) spin state. Therefore, 
neither of these spin states participate in the total spin 
channel photoassociation if they coexist in the con- 
densate only with an m/ = spin state. To illustrate 
these predictions, PA spectra across the lines e and are 
taken with condensates containing different mixtures of 
spin states as shown in Fig. [2] 

The data shown in Figs. [Ha) and dfb) are taken with 
similar conditions to Fig.Q] whereas in Fig.[2][c) and[2ld) 
the photoassociation power is decreased to 1.5 mW to 
reduce the mechanical effects on the m/ = — 1 atoms 
caused by the rapid loss of the m/ = atoms from the 
trap. These spectra are taken with frequency steps of 2 
MHz. The observed data in this figure shows that the 
rrif = — 1 and mf = 1 spin states participate in the pho- 
toassociation for the lines e and if both of them coexist 
in the condensate, whereas if m/ = —1 spin state coex- 
ists with the mf — spin state it does not participate in 
the photoassociation process. This observation matches 
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FIG. 2: Photoassociation spectroscopy of e and # with dif- 
ferent mixtures of spin states, (a) and (b) Lines e and 6 for 
mixture of to/ = —1 and to/ = 1 spin states, (c) and (d) for 
mixture of to/ = — 1 and to/ = 0. The observed data points 
corresponding to m/ = ( — 1,0, 1) spin states are represented 
by circles, triangles, and squares respectively. 



FIG. 3: Photoassociation spectroscopy through lines a, r), fi, 
and £ for an to/ = spin state condensate to an F — 2 molec- 
ular state. The data points are taken with the polarization 
vector of the PA light along and perpendicular to an external 
magnetic field (squares) and (circles). 



our analysis based on scattering channels and confirms 
that F = 1 for lines e and 9. 

To identify F = 2 PA lines we note that a AF = 0, 
AMp = transition from a total two atom collisional 
spin state |2, 0) to a bound molecular state |2,0) is for- 
bidden by the electric dipole selection rule. Therefore 
the lines with F = 2 are suppressed for an mt =0 spin 
state condensate with 7r-polarized PA light. Such lines 
are found by taking a PA spectrum of a pure m/ = 
condensate in the presence of an external B-field of 1 G 
along the polarization of the PA light. The observed data 
for this case, shown in Fig.[3j indicates that the lines a, 
77, y,, and £ are absent. By changing the polarization of 
the PA light perpendicular to the quantization axis de- 
fined by the applied B-field, all of these lines reappear in 
the spectrum. A comparison between the observed data 
with polarization aligned parallel (square) and perpen- 
dicular (circle) to the external B-field is shown in Fig.[3l 
which confirms F = 2 as a good quantum number for 
these lines. Fig.[3£c) also indicates that the line fj, is split 
into two F = 2 components. 

From these observations and the spacing of the lines we 
are able to assign F and / for the observed lines. This is 
accomplished by noting that for possible values of F = 1 , 
2, and 3 the eigenvalues of F 2 are 2, 6, and 12. Therefore 
the separation between molecular states with F = 1 and 
2 should be 2 /3 of the frequency separation between the 
lines with F = 2 and 3. F numbers for the remaining 
lines are readily deduced by their spacing and order in 
each hyperfine grouping. The total nuclear spin for these 
lines is found by decomposing the initial scattering wave 
function of the two colliding atoms. This analysis shows 
that the total nuclear spin must be 1 for the total spin 
scattering channel and predominantly / = 3 for the 
total spin 2 scattering channel. TableQ] gives the assigned 
values of F and / for the observed lines. 



TABLE I: Assigned total angular momentum and nuclear spin 
of the molecular states appearing in the PA spectra of Fig.Q] 
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For the l g potential, the molecular hyperfine constant 
is comparable to the small rotational constant of the 
state. Therefore, the effective Hamiltonian for the l g 
rotation-hyperfine structure can be written as [iij ]. 



Hh-r — a v Ltl 

= a v I g n 
+ B V {F 2 ^ 



B V J 



P-2F- Z I- z 



F+f_-F_f+}, (1) 



where the operators are defined with respect to the inter- 
nuclear axis z. Since changing the signs of F s , Is, and 
f2 results in a degenerate state, the states are labeled by 
F, /, I, and i, where f — F s ■ Q and i = lg ■ fl. The 
eigenvalues of the Hamiltonian are obtained by diagonal- 
izing its corresponding matrix in the \F,Fg,I,Ig) basis 
set. The coupling is sufficiently weak that labeling the 
final states by the almost good quantum numbers of / 
and i is justified. These eigenvalues are fit to the ob- 
served spectrum of the m/ = —1 condensate of Fig.Q] 
using the least squares approach with a v , B v , and the 
frequency offset w ffsct as fitting parameters. In case of 
the split lines, the stronger one is used for fitting. Us- 
ing the aF=2 lines the fitting specifies a v and B v to be 
667(5) and 27.6(1.3) MHz, respectively. The positions of 
the eigenvalues using these parameters are shown as a 
stick spectrum in Fig. [T] which is in good agreement with 
the observed locations of the atom loss lines in the PA 
spectrum. The stick spectrum is labeled by the domi- 
nant part of the corresponding eigenvectors which is in 
complete agreement with the assigned quantum numbers 
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FIG. 4: Calculated values for Aa and for the PA limited 
lifetime of BEC. The dotted line is the nominal value of Aa. 
The intensity here is 20 times that used to take the spectra 
in Fig.[H 



in Tabled] 

A few general comments should be made for the ob- 
served as well as predicted spectra. First, the molecular 
state corresponding to the |2, —1, 3, —2) eigenstate is not 
observed under any experimental condition. This line is 
also absent in a scan of the v — 153 vibrational level. 
Second, predicted locations of the lines e and 9 corre- 
sponding to the aF=o scattering channel are in poorer 
agreement with their observed locations compared to the 
other lines. These two lines are connected by arrows to 
their assumed eigenvalues in Fig.[TJ There appears to be a 
discrepancy both in the separation of the lines and their 
apparent origin for the hyperfine-rotation Hamiltonian. 
It was unclear as to whether different fit parameters or a 
revision to the Hamiltonian is needed to correct the de- 
viations. Last, we note that some of the lines appear to 
be split on the order of 30 to 160 MHz. These split lines 
are connected to their labels in Fig.Q] for clarity. The 
Hamiltonian presented in Eq. (JTJ) predicts that all lines 
should be doubly degenerate and does not account for 
this splitting. This points to an additional interaction 
not considered in this simple model. These discrepancies 
are important for refinement of the molecular potential 
theory. 

One of our motivations for studying this system is to 
assess the suitability of photoassociation induced change 
in two s-wave scattering lengths, Aa = ap = 2 — of=o, in 
order to manipulate spin-dependent properties of a con- 
densate. To determine the variation of Aa, the observed 
data of Fig. [1] is fitted to a theoretical formula for the in- 
elastic loss rate (K- lne i) to find the width and amplitude 
for each line as discussed in Ref. [ll|- These fit param- 
eters along with the condensate density (no) are then 
used to calculate the estimated change in Aa as shown 
in Fig.HJa). 

We also calculate a photoassociation limited lifetime, 
tpa = (-Kmci x no)~ 1 , during which spin dependent dy- 
namics can be observed before the condensate is depleted 



by the PA atom loss as shown in Fig.[4jb). The plots in 
Fig. [4] assume a PA intensity 20 times higher than that 
used to measure the spectrum. One of the more inter- 
esting frequency ranges in Fig. [4] is between lines 5 and 
e which are associated with different scattering channels. 
In this frequency range the influences of the two lines 
reinforce each other for a large change in Aa, but with 
only a modest reduction in the 3 s lifetime of the con- 
densate due to the PA light. From this data it is clear 
that the value of Aa can be altered significantly from 
its nominal value of 1.45(32)as 7] up to a total cance- 
lation, with sufficient PA limited lifetimes needed to see 
the change in spinor dynamics. However to change the 
sign and hence the magnetic nature of the condensate, 
the values needed correspond to challengingly short life 
times compared to the time needed to relax to the mag- 
netic ground state. Hopefully, other molecular states will 
have a slightly more favorable width, amplitude, or spac- 
ing. 

In conclusion, we performed photoassociation spec- 
troscopy on a spin-1 condensate. The spin dependent 
photoassociation spectrum is used to identify good quan- 
tum numbers for some of the molecular states, which are 
in agreement with theoretical predictions of a hyperfine- 
rotation Hamiltonian. The spectrum is also used to pre- 
dict the photoassociation limited lifetime and change in 
Aa. It is shown that optical Feshbach resonances to 
a molecular state with hyperfine structure are a viable 
method of altering the spin-dependent mean-field interac- 
tion energy to change the spinor dynamics of the conden- 
sate. The results of this study provide a valuable test case 
to answer general questions about how to model molec- 
ular potentials in the presence of hyperfine and rotation 
interactions. A high resolution photoassociation study of 
other vibrational levels will provide a deeper understand- 
ing of molecular potentials. This work is supported by 
the National Science Foundation, PHYS-0605049. 
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